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The NAVSTAR Global Positioning System is in- 
tended to provide observers with instantaneous 
measurements of position and velocity in three di- 
mensions to 10 m accuracy in 1985- In order to 
obtain position in three dimensions, "pseudo-range" 
measurements are computed to four satellites on the 
basis of the difference in the time of transmission 
and the time of receipt of signals transmitted by 
the satellite. Since the time of receipt of the 
signals is biased due to the error in the observer's 
clock, the computed ranges are called "pseudo-ranges 
and the measurements are made to four satellites so 
that the error in the observer's clock as well as 
the three components of position can be computed. 

In order to test the concept in 1979 j six satellites 
are being launched into orbits which will allow an 
observer in the Southeastern part of the United 
States to receive simultaneous data from four satel- 
lites for several hours each day. Both test and 
operational satellites will be in circular twelve 
hour orbits at an inclination of 63 degrees. The 
operational system will have eight satellites in 
each of three orbit planes to provide world-wide 
continuous coverage. 

A number of proposals have been made to utilize 
the GPS satellites to obtain positions to a few 
centimeters accuracy for use in crustal motion re- 
search. In each instance, orbit and satellite clock 
errors are minimized by computing the relative posi- 
tion of ground sites. Because of the high altitude 
of the satellites, these errors will not be signi- 
ficant even for stations separated by hundreds and 
perhaps thousands of kilometers, depending on 
whether predicted, fitted, or specially computed 
ephemerides are used (Anderle, in press). MacDoran 
(in press) proposes to make measurements on the 
spread spectrum signal broadcast by the satellite. 
Shapiro (private communication) would make measure- 
ments on an additional signal source he proposes 
to place aboard the satellites. Measurements 
could also be made on side-tone signals if these 
were provided by the satellite. In this report I 
will address use of the signals transmitted by 
the satellite and recorded by prototype equip- 
ment which' is now being tested. Each of these 
approaches has certain advantages and disad- 
vantages : 
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The satellite transmits on two L band frequencies 
centered at 1575 (Lq) and 1227 (lg) Mhz in order to 
provide a capability for correction for first order 
ionospheric effects. The carriers are fully sup- 
pressed by a bi-phase pseudo-random modulation at 
a 10 megabit rate. This "P code" is supplemented by 
a lower bit rate, shorter code sequence, "C/A code" 
on the Lj signal which permits rapid acquisition 
of the signal. The code is also generated by the 
receiver, where it is shifted in time to obtain 
maximum correlation with the received signal. The 
time of maximum correlation is recorded as the time 
of receipt of the signal. Subtraction from broad- 
cast time of transmission of the signal gives the 
apparant travel time to the signal, and multiplica- 
tion of the difference by the velocity of light 
yields the pseudo range to the satellite. Once 
the correlation of the codes is made, the carrier 
signal is reconstructed and Doppler measurements 
axe made. 

Eight receiver channels would be required to 
maintain continuous lock on the signal from 
four satellites as well as the (i^-I^) signals. 

To reduce receiver size, power, weight and cost, 
fewer channels are used. Since the initial search 
for the signal correlation can take tens of seconds, 
five channels is the minimum number used when high 
data rates are desired; four of the channels main- 
tain continuous lock on the signals from four 
satellites while the fifth channel rapidly se- 
quences through the (Lj-Lg) signals. Simpler 
receivers have two channels which lock on the 
L-^ and Lg signals from a sateliite, then switch 
to the next satellite, or a single channel which 
sequences through the Lj then Lg signals on each 
satellite in turn. The lowest cost receiver uses 
a single channel which sequences through only the 
Lq C/A signal on each satellite. The primary 
output of these receivers, depicted in figure 1, 
is the range obtained by correlation of the 
pseudo random noise signal. Doppler is counted 
to obtain velocity and reduce the noise of the 
PRN measurements. The NAVSTAR Geodetic Receiver 
(NGR) built by the Naval Surface Weapons Center 
and Stanford Telecommunications Incorporated 
(figure 2) has two channels and records range at 
Lq and Doppler at Lj and L2, sequencing through 
satellites according to computer commands. 

The measurement precision depends on antenna 
gain and the receiver hardware. Examples of pre- 
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FIGURE 1 


oision for two antenna gains are given below 
based on the equations given in figure 3 for 
monitor stations (similar to X sets) used to ob- 
tain data for orbit computations; the precision 
is also given for the NGR. 





Table 1 


Sample Measurement Precisions 


Equipment 

Monitor Station 

NGR 


(6 sec . average ) 


Range Precision 

0 db antenna 

36 cm 

65 cm 

21 db antenna 

0.1* cm 

5 cm 

Doppler Precision 

0 db antenna 

O.ll* cm 

0.17 cm 

21 db antenna 

0.08 cm 

0.01 cm 


The above figures are based upon a receiver power 
of 163 dbw, which is 3 dbw poorer than that ob- 
served for the first two satellites launched. 

The zero db antenna is a small stub which does 
not have to be aimed. A 21 db antenna is equi- 
valent to a one meter dish which must track each 
satellite, although the required aiming accuracy 
is only about 10 degrees. If the measurement 
precision at each of the L— band frequencies is 
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the same, the precision of the range corrected for 
ionospheric refraction will be a factor of four 
worse than that fundamental precision after ap- 
plication of the correction algorithm. However, 
by tracking the (L^-Lj) the refraction correction 
can be obtained without significant degradation 
of the basic precision, as shown in figure 3. 

The precision inferred from monitor station data 
which has geen aggregated to 15 minute intervals 
is currently 20 cm. This is higher than ex- 
pected for aggregated data, according to table 1; 
however, the value of 20 cm was obtained as the 
residuals of fit to five to ten hours of data. 
Figure 5 shows that oscillator variations in- 
tegrate to range errors of 50 to 100 cm after 
5 to 10 hours, which easily accounts for the re- 
sult obtained.' Hermann (private communication) 
analyzed 15 minute segments of data obtained on 
two satellites during an initial demonstration 
of the NGH. He found the noise level of the 
range measurements matched the expected 65 cm 
for one satellite. The noise of frequency 
measurements on each satellite at matched 
that expected for the frequency counter used in 


the demonstration, although the equipment was 
only capable of validating the precision of the 
phase measurements to the 20 cm level. Fre- 
quency measurements at Lp were considerably 
worse than expected, probably due to interf er- 
ring signals encountered in the area. 

The difference in positions of two stations 
based on pseudo-range measurements to four satel- 
lites is equivalent in information content to 
an inter -ferometric determination of relative 
station positions. In such solutions, the pre- 
cision of the position determination is a factor 
of two to three- worse than the precision of the 
range measurements due to the Geometric Dilution 
of Precision (GD0P) arising from the simultaneous 
solution for position and time bias. Thus the 
expected precision for solutions for relative 
station positions based on X-set dat a ove r a 
six hour time span would be J" 2x3x35/73^00 =2.5 cm, 
where /2 arises from differencing data from two 
stations, 3 is a GD0P, 36 cm is the precision of 
six second data, and 3600 is the number of obser- 
vations in a six hour interval. About the same 
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where 

B c = code channel bandwidth - 25 hz 

S/W Q = signal to noise ratio in hz - lO^^^hz 


c = velocity of light 

f. = code frequency -10.23 mhz 

Tp= L-^ code loop time constant - 10 sec 

M = number of operating channels 

T v = nominal sampling interval ~ 10 sec 

q = c/64 f = l/64th of a code chip 

T 12 = m i n ( T jl2 > T Ts/3^’ where T 2 = l 1 -Ij 2 code loop constant ~ 160 sec, T Tg = time slot of MS 
tracking sequence - 60 sec 
Bl = costas loop bandwidth for Doppler -3.33 hz 

f c = carrier frequency -1575*42 mhz 


q c = c/61* f c = l/64th of a carrier cycle 
t = integrated doppler measurement interval ~ 6 sec 

Tj = time interval per unit change in range due to ionosphere '197 sec/m 
interface control document 12436, MCS/NSWC 


FIGURE 3 
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level of accuracy was obtained by Goad (private 
communication) in a computer simulation which 
considered the geometric and other factors in 
the solution in a more realistic fashion. 

Because of the higher precision of phase 
measurements and the successful use of such 
measurements of Navy Navigation Satellites in 
geodetic solutions, the Naval Surface Weapons 
Center is exploring the use of such data from 
GPS satellites for positioning. In this approach 
phase measurements and integral Doppler counts 
are made over successive time intervals during 
the passage of a satellite between the start 
and end of the interval can be computed. For 
the lower altitude, polar. Navy satellites, 
the range difference observations obtained dur- 
ing one pass of the satellite over the station 
provide a determination of station latitude and 
the range to the satellite at the time of closest 
approach of the satellite to the station. Data 
from two passes at different longitudes are 
required to resolve the two slant range measure- 
ments into station longitude and height to 
acceptable precision. For the higher altitude 
GPS satellites, a three-dimensional solution 
for station position can be obtained from obser- 
vations made on a single pass of the satellite 
over the station, but not to the precision pos- 
sible using data from two passes. The range 
differences observed during a pass can be 
treated in two ways: the range differences 

over successive time intervals can be assumed 


to be uncorrelated, or, if the cycle count is 
not reset at the time of readout, the range 
differences can be accumulated to provide a 
series of range measurements subject to a 
range bias common to the observation set. 
Simulations by Anderle (in press) indicate that 
biased range data during a pass yields a pre- 
cision in the determination of two components 
of station position which is about equal to the 
precision of the range difference measurement. 
Uncorrelated range difference data provides a 
precision which is a factor of three or more 
worse than the precision of measurement. Due 
to the high measurement precision possible, 
this would still be a highly useful data type 
which would permit observations of two or more 
satellites sequentially with a two channel re- 
ceiver. However, range difference data, 
treated as either correlated or uncorrelated, 
set severe requirements on the oscillator used 
to make the measurements. Figure 1 (Fell, 
private communication) shows that a Cesium 
oscillator will produce systematic errors which 
would reach 20 cm during a fit of mean frequency 
offset to eight hours of range difference data. 
The figure is a sample result for one sequence 
of data. Fell also gives the root mean square of 
the error for a large set of sequences of such 
data in figure 5 as a function of time within the 
pass; the rms error ranges frcsn 10 to 15 cm, with 
the maximum error in the center and at the ends 
of the pass. The effect of this systematic error 
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where 

c = velocity of light 

Bj-(DLL) = one-sided loop noise bandwidth = 2 hz 

g (PLL) = one-sided phase loop bandwidth = l6 hz 

S/N q = signal to noise ratio in hz (10^ B ^ BB ) 

f = code frequency~10.23 mhz 

Bj = IF, bandwidth = 10^ hz 

f = carrier frequency - 157 ^- ^2 mhz 
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FIGURE It 


56 




FIGURE 5 


on station position is being computed. But it 
will clearly dominate the effects of the measure- 
ment error. 

To summarise, pseudo-range measurements to 
four GPS satellites based on correlation of the 
pseudo random code transmissions from the satel- 
lites can be used to determine the relative po- 
sition of ground stations which are separated 
by several hundred kilometers to a precision at 
the centimeter level. This precision is attained 
within 12 hours using a small antenna or much 
more rapidly with a directional antenna. Carrier 
signal measurements during the course of passage 
of two satellites over a pair of stations would 
also yield centimeter precision in the relative 
position, but oscillator instabilities would 
limit the accuracy by an as yet undetermined 
amount. Measurement precisions of code and 
carrier signals have generally been consistent 
with test conditions, but have not been tested 
to design levels yet. The accuracy of solutions 
based on either type of data would be limited by 
unmodeled tropospheric refraction effects which 
would reach five centimeters at low elevation 
angles for widely separated stations. 
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